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A  thermodynamic  study  of  Glycerol  Steam  Reforming  (GSR)  for  hydrogen  production  with  in  situ  carbon 
dioxide  and  hydrogen  (reaction  products)  simultaneous  removal  was  performed.  The  sorption-enhanced 
membrane  reactor  (SEMR)  was  divided  into  multiple  sub-Gibbs  reactors  and  the  Gibbs  free  energy 
minimization  method  was  employed.  The  effects  of  temperature  (600-800  K),  molar  water-to-glycerol 
feed  ratio  (WGFR)  (3-9),  pressure  (1-5  atm)  and  fraction  of  hydrogen  and  carbon  dioxide  removal  (/, 
0-0.99)  on  the  GSR  process  were  target  of  investigation.  A  hydrogen  yield  (total  moles  of  hydrogen 
produced/mole  of  reacted  glycerol)  very  close  to  the  stoichiometric  value  of  7  was  obtained  at  700  K, 
WGFR  of  9, 1  atm  and  for  /Cq2  =  0.99  and  /Hz  =  0.80.  This  corresponds  to  an  enhancement  of  217%,  47% 
and  22%  in  terms  of  hydrogen  yield  comparatively  to  the  traditional  reactor  (TR),  sorption-enhanced 
reactor  (SER)  with  carbon  dioxide  capture  (fCo2  =  0.99)  and  membrane  reactor  (MR)  with  hydrogen 
separation  (fHl  =  0.80),  respectively.  In  terms  of  coke,  its  formation  was  only  observed  under  WGFRs 
below  the  stoichiometric  value  of  3. 


©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  use  of  glycerol,  the  main  by-product  of  the  biodiesel 
manufacture  process,  for  hydrogen  production  is  not  only  envi¬ 
ronmentally  more  attractive  than  fossil-based  routes,  but  also 
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valorizes  glycerol  itself  thus  making  biodiesel  more  competitive. 
The  Glycerol  Steam  Reforming  (GSR)  is  an  endothermic  and 
equilibrium-limited  process  that  requires  high  operating  temper¬ 
atures  (Eq.  (1)  -  cf.  Table  1).  The  GSR  process  involves  the 
decomposition  of  glycerol  (Eq.  (2)  -  cf.  Table  1)  followed  by  the 
water-gas  shift  (WGS,  Eq.  (3)  -  cf.  Table  1 )  reaction  (multiplied  by  a 
factor  of  3). 

Even  though  a  theoretical  production  of  7  mol  of  hydrogen  per 
mole  of  consumed  glycerol  should  be  observed,  the  existence  of 
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Notation  and  glossary 

List  of  acronyms 

GSR 

Glycerol  Steam  Reforming 

MR 

membrane  reactor 

RWGS 

reverse  water-gas-shift 

SER 

sorption-enhanced  reactor 

SEMR 

sorption-enhanced  membrane  reactor 

SRM 

steam  reforming  of  methane 

TR 

traditional  reactor 

WGFR 

water-to-glycerol  feed  ratio 

WGS 

water-gas  shift 

side  reactions  (e.g.  methanation)  leads  to  a  decrease  in  hydrogen 
production.  Theoretical  studies  on  new  reactor  configurations  that 
combine  GSR  and  carbon  dioxide  [1,2]  or  hydrogen  3]  selective 
removal  (reaction  products  in  Eq.  (1))  have  allowed  to  conclude 
that  such  intensified  processes  permit  to  enhance  the  hydrogen 
production  and  decrease  the  production  of  both  methane  and 
carbon  monoxide  by-products  by  shifting  the  thermodynamic 
equilibrium.  Moreover,  such  processes  allow  the  attainment  of 
similar  or  even  better  performance  than  in  a  traditional  reactor  (TR) 
while  operating  at  milder  operating  conditions.  However,  the  effect 
of  the  selective  removal  of  both  products  simultaneously  during 
GSR  has  never  been  studied.  Such  a  process  could  be  achieved  with 
a  multifunctional  Sorption  Enhanced  Membrane  Reactor  (SEMR), 
which  would  simultaneously  carry  out  the  GSR  reaction  and 
remove  pure  hydrogen  and  carbon  dioxide  from  the  reaction  zone 
by  means  of  a  hydrogen  selective  membrane  and  a  carbon  dioxide 
selective  sorbent,  respectively  (Fig.  1).  The  continuous  use  of  such 
reactor  configuration  would  require  two  parallel  reactors,  being 
that  while  one  of  them  is  producing  pure  H2  through  GSR  (which 
exits  in  the  permeate  stream  and  is  not  mixed  with  any  other 
species,  apart  from  the  inert  carrier  gas,  due  to  membrane  selective 
permeation  -  left  side  of  Fig.  1 ),  the  other  is  being  regenerated,  i.e., 
carbon  dioxide  is  being  removed  and  exits  in  the  retentate  stream 
(again  N2  can  be  used  as  purge  gas  -  right  side  of  Fig.  1).  In  other 


Table  1 

Reactions  in  equilibrium  considered  during  the  GSR  thermodynamic  simulations. 


Reaction 

AH;?98  K  (kj  mol  :) 

Reaction  number 

C3H803  +  3H20  ^  7H2  +  3C02 
(GSR  —  Glycerol  Steam  Reforming) 

128 

(1) 

C3H803  ^  3 CO  +  4H2 
(decomposition  of  glycerol) 

251 

(2) 

CO  +  H20  ^  h2  +  C02 
(WGS  —  water-gas  shift) 

-41 

(3) 

CO  +  3H2  —  CH4  +  H20 
(methanation) 

-206 

(4) 

C02  +  4H2  CH4  +  2H2Oa 

-165 

(5) 

2CO  +  2H2  CH4  +  C02b 

(dry  reforming  of  methane) 

-247 

(6) 

2CO  C02  +  C 

(Boudouard  reaction) 

-172 

(7) 

CH4  ^  2H2  +  C 
(methane  cracking) 

75 

(8) 

CO  +  H2  ^  H20  +  C 
(carbon  monoxide  reduction) 

-131 

(9) 

C02  +  2H2  ^  2H20  +  C 
(carbon  dioxide  reduction) 

-90 

(10) 

a  Sum  of  both  reverse  of  Eq.  (3)  and  Eq.  (4). 
b  Sum  of  Eqs.  (3)  and  (4). 


words,  each  reactor  has  two  operation  modes  being  that  along  time 
each  of  them  goes  through  the  reaction-regeneration  cycle 
repeatedly.  Since  the  goal  is  to  shift  the  thermodynamic  equilib¬ 
rium  of  GSR  during  reaction  mode  by  retaining  C02  in  the  sorbent, 
once  the  sorbent  gets  saturated  and  C02  starts  breaking  through 
the  column  the  reaction  is  ended  (Fig.  1(a)).  At  this  point,  the  set  of 
valves  moves  this  reactor  feed  (glycerol  and  steam)  to  e.g.  N2  so  that 
the  sorbent  bed  can  be  regenerated  (Fig.  1(b)).  During  this  last  step 
no  H2  is  produced  in  this  reactor,  being  produced  in  the  other 
reactor  that  is  now  in  the  reaction  stage.  Moreover,  although  C02 
might  affect  membrane  permeability  towards  H2,  this  is  not  very 
pronounced  as,  for  instance,  it  occurs  with  CO,  and  such  effect  is 
reversible  (thus,  as  long  as  C02  is  removed,  high  permeability  to¬ 
wards  H2  is  restored)  [4].  The  main  requirement  is  that  operating 
conditions,  namely  in  terms  of  temperature,  are  compatible  for  the 
GSR  catalyst,  C02  sorbent  and  H2-selective  membrane,  which  is  the 
case  [5  . 

In  this  work  an  equilibrium  thermodynamic  analysis  of  GSR  in  a 
SEMR  was  done,  for  the  first  time,  with  ASPEN  PLUS  software  and 
the  effect  of  the  operating  temperature,  pressure,  water-to-glycerol 
feed  ratio  (WGFR)  and  H2  removal  fraction  (for  constant  C02 
removal)  on  hydrogen  production  was  accessed.  Moreover,  a 
comparison  between  the  results  obtained  under  equilibrium  con¬ 
ditions  for  the  SEMR  and  for  the  other  types  of  reactors  (TR, 
membrane  reactor  (MR)  and  sorption  enhanced  reactor  (SER))  was 
done  as  well.  Finally,  the  limits  in  terms  of  WGFR  for  avoiding  coke 
formation  at  different  temperatures  was  discussed.  This  thermo¬ 
dynamic  analysis  is  crucial  as  it  may  indicate  the  suitability  of 
applying  such  intensified  reaction-separation  process  for  hydrogen 
production  from  glycerol.  Moreover,  and  although  thermodynamic 
considerations  alone  do  not  provide  evidence  for  the  practicability 
of  the  design,  it  provides  valuable  information  about  the  limits  that 
can  be  attained  in  actual  operation,  and  gives  indication  on  how  to 
proceed  and  which  conditions  have  to  be  applied  in  order  to  ach¬ 
ieve  maximum  performance. 

2.  Methodology 

ASPEN  PLUS  V7.3  was  used  to  perform  all  the  simulations.  The 
Gibbs  free  energy  minimization  methodology  (nonstoichiometric 
method),  which  calculates  the  equilibrium  composition  through 
the  minimization  of  the  Gibbs  free  energy  of  a  specific  set  of  species 
was  applied  (a  more  detailed  description  of  the  methodology  is 
available  in  Supplementary  data). 

For  such  simulations,  the  species  included  were  hydrogen,  car¬ 
bon  dioxide,  carbon  monoxide,  methane,  glycerol,  water  and  solid 
carbon.  By  considering  these  species,  some  of  the  possible  reactions 
that  might  be  in  equilibrium  together  with  GSR  (Eqs.  (2)  and  (3)), 
are  the  ones  shown  in  Table  1.  The  reaction  in  Eq.  (5),  which  is  the 
sum  of  both  reverse  of  Eq.  (3)  and  Eq.  (4),  is  included  in  Table  1  with 
the  single  purpose  of  facilitating  the  analysis  of  the  stoichiometric 
relations  in  the  SEMR. 

The  modular  approach  used  to  simulate  the  SEMR 
(Supplementary  data,  Fig.  SI)  is  represented  by  (n  +  1)  sub¬ 
reformers  and  n  sub-separators.  Each  sub-reformer  consists  on  a 
Gibbs  reactor  which  performs  the  thermodynamic  equilibrium 
calculations  through  the  method  previously  explained.  The  sub¬ 
separators  consist  on  component  separators  that  separate  chemi¬ 
cal  species  based  on  specified  flows  or  split  fractions,  which  in 
practice  is  related  with  factors  like  membrane  and  sorbent  selec- 
tivities,  membrane  area  and  sorbent  quantity,  etc.  For  n  =  0  a 
simple  Gibbs  reactor  is  simulated;  for  n  =  1  a  single  separator  and 
two  sub-reformers  are  used;  and  so  on.  The  value  of  n  depends  on 
the  values  of  the  global  removal  fractions  of  C02  (fco2)  and  H2  (fHz) 
specified  (higher  values  of/ imply  higher  values  of  n).  It  is  worth 
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Fig.  1.  Schematic  view  of  the  conceived  SEMR  based  on  2  parallel  reactors  configuration  for  continuous  operation  and  corresponding  outlet  concentrations  histories  in  the  retentate 
stream  during  (a)  reaction  and  (b)  regeneration  stages. 


mentioning  that  a  similar  methodology  has  been  used  in  previous 
works  [1,3]. 

The  thermodynamic  analysis  was  performed  for  temperatures 
in  the  range  of  600-800  K,  pressures  between  1  and  5  atm  and 
WGFRs  between  3  and  9  (realistic  operation  conditions  for 
hydrogen  selective  membranes  [6,7]  and  carbon  dioxide  sorbents 
[5,8,9  ,  which  are  close  to  those  typically  employed  in  the  GSR 
[10,11  );  the  study  considered  also  hydrogen  removal  fractions  in 
the  range  of  0-0.99  and  a  carbon  dioxide  removal  fraction  of  0.99 
(for  the  reasons  explained  in  a  posterior  section).  For  all  the  cases, 
the  yields  of  hydrogen,  methane,  carbon  dioxide  and  carbon 
monoxide  were  analyzed.  These  were  defined  as  follows: 

=  (jl[,  separated  +  ^i,  unseparated)  /^converted  glycerol  (H) 

where  Y[  is  the  yield  of  species  i  and  nconverted  glycerol  is  the  molar 
flow  rate  of  glycerol  that  is  converted  into  all  reaction  products. 

3.  Results  and  discussion 

3.2.  GSR  with  in  situ  hydrogen  or  carbon  dioxide  removal  - 
comparison 

Before  performing  the  thermodynamic  equilibrium  calculations 
for  the  SEMR,  the  thermodynamics  of  GSR  performed  in  a  TR,  SER  or 
MR  was  analyzed.  In  order  to  fully  understand  the  effect  of  the 
simultaneous  removal  of  both  hydrogen  and  carbon  dioxide  during 
GSR  on  the  thermodynamic  equilibrium,  it  is  essential  to  analyze 
first  the  effect  of  the  individual  removal  of  either  hydrogen  or 
carbon  dioxide  on  the  GSR  thermodynamics.  A  comparison  be¬ 
tween  both  individual  effects  (hydrogen  and  carbon  dioxide 
removal),  which  to  the  best  of  our  knowledge  has  never  been  done 
before,  is  included  as  well. 

Complete  conversion  of  glycerol  was  observed  for  all  the  range 
of  temperatures  analyzed  (data  not  shown),  in  agreement  to  what 
was  observed  by  other  authors  1  — 3,12,13 ].  This  means  that  in  the 
subsequent  analysis  Eq.  (2)  is  almost  not  taken  into  account  as  it  is 
completely  shifted  towards  reaction  products.  Coke  formation  was 


not  observed  for  all  the  conditions  used  in  this  section,  so  Eqs. 
(7)-(10)  will  now  be  discarded. 

3.2.2.  Traditional  reactor  (ft  =  0) 

The  variation  with  temperature  of  the  dimensionless  (normal¬ 
ized  with  the  molar  flow  rate  of  glycerol  fed  to  the  reactor)  molar 
flow  rate  of  glycerol,  water,  hydrogen,  carbon  dioxide,  carbon 
monoxide  and  methane  in  the  thermodynamic  equilibrium  of  a  TR 
is  presented  in  Fig.  2. 

In  terms  of  hydrogen  and  carbon  dioxide  dimensionless  molar 
flow  rates,  it  can  be  seen  that  they  increase  with  temperature  up  to 
a  certain  point,  stabilize  and  then  suffer  a  slight  decrease  (Fig.  2).  In 
order  to  understand  this  behavior,  the  reactions  in  Eqs.  (3)  and  (4) 
should  be  considered.  At  lower  temperatures  (-600  K),  the  strongly 


Fig.  2.  Effect  of  temperature  on  the  thermodynamic  equilibrium  composition  of  a 
traditional  reactor  at  WGFR  of  9  and  1  atm. 
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exothermic  methanation  reaction  (Eq.  (4))  is  favored,  in  agreement 
with  the  maximum  methane  and  water  and  minimum  carbon 
monoxide  dimensionless  molar  flow  rates  observed.  Particularly, 
the  dimensionless  molar  flow  rate  of  water  for  such  temperature  is 
slightly  higher  than  the  value  of  9  in  the  feed,  thus  indicating  water 
formation.  When  temperature  increases  the  methanation  reaction 
is  more  disfavored  (progressive  inhibition  of  methane  production 
as  compared  to  WGS  reaction  (less  exothermic  -  Eq.  (3))),  which 
progressively  becomes  dominant  (CO2  and  H2  production  and  H2O 
consumption).  The  stabilization  of  the  hydrogen  and  carbon  diox¬ 
ide  production  at  higher  temperatures  is  related  to  the  highly 
enhanced  reverse  WGS  (RWGS),  which  leads  to  an  increase  in  the 
water  production.  The  faster  stabilization  of  the  carbon  dioxide 
production  observed  (maximum  at  lower  temperature  than 
hydrogen)  is  due  to  the  fact  that  while  hydrogen  is  being  consumed 
through  RWGS  and  saved  through  the  inhibition  of  methanation 
(3  mol  of  H2  per  mole  of  CH4  unproduced),  carbon  dioxide  is  only 
being  consumed  through  RWGS.  When  the  temperature  is  high 
enough  so  that  there  is  no  methane  formation,  the  yield  of 
hydrogen  starts  decreasing  and  the  RWGS  becomes  the  only 
dominant  reaction  (hydrogen  and  carbon  dioxide  decrease  at  a 
similar  rate).  Since  carbon  monoxide  is  consumed  through  WGS 
and  methanation,  which  are  both  exothermic  reactions,  higher 
temperatures  lead  to  higher  carbon  monoxide  production.  For  the 
conditions  here  analyzed  the  material  balances  performed  for  each 
species  considering  only  stoichiometry  of  reactions  2,  3  and  4  were 
all  consistent  with  the  simulations  results,  thus  indicating  that 


there  are  no  other  reactions  occurring  besides  the  ones  previously 
indicated. 

Regarding  the  dry  reforming  of  methane  (Eq.  (6)),  suggested  by 
Adhikari  et  al.  [12]  to  be  responsible  for  the  behavior  of  carbon 
dioxide  production,  once  it  depends  on  both  WGS  and  methanation 
(Eqs.  (3)  and  (4)),  which  have  been  target  of  discussion  already,  it  is 
implicitly  already  analyzed  here. 

3.2.2.  Sorption-enhanced  reactor  versus  membrane  reactor 

The  variation  of  the  yield  of  hydrogen,  methane,  carbon  mon¬ 
oxide  and  carbon  dioxide  with  temperature  and  carbon  dioxide 
removal  fraction  (in  the  SER)  and  hydrogen  removal  fraction  (in  the 
MR)  is  presented  in  Figs.  3  and  4,  respectively.  Besides  increasing 
with  temperature,  the  yield  of  hydrogen  increases  with  the  fraction 
of  carbon  dioxide  or  hydrogen  removal  as  well.  It  can  be  observed 
by  comparing  Figs.  3(a)  and  4(a)  that  higher  hydrogen  production 
can  be  attained  by  selectively  removing  hydrogen  from  the  reaction 
medium  instead  of  carbon  dioxide,  especially  at  lower  tempera¬ 
tures.  This  difference  is  associated  to  the  inhibition  of  the  formation 
of  the  main  GSR  by-product,  methane.  When  carbon  dioxide  is 
removed  from  the  reactor,  the  equilibrium  of  the  WGS  reaction  (Eq. 
(3))  is  shifted  towards  the  production  of  more  carbon  dioxide,  and 
inherently  hydrogen,  and  carbon  monoxide  consumption.  Since 
carbon  dioxide  does  not  participate  in  the  reversible  methanation 
reaction  (Eq.  (4)),  its  removal  does  not  directly  affect  this  reaction. 
Instead,  it  is  only  affected  by  the  enhanced  production  of  hydrogen 
and  carbon  monoxide  consumption  (opposite  effects).  On  the  other 
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Fig.  3.  Effect  of  temperature  and  carbon  dioxide  removal  fraction  of  the  SER  on  the  yield  of  (a)  hydrogen,  (b)  methane,  (c)  carbon  monoxide  and  (d)  carbon  dioxide,  at  WGFR  of  9 
and  1  atm. 
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Fig.  4.  Effect  of  temperature  and  hydrogen  removal  fraction  of  the  MR  on  the  yield  of  (a)  hydrogen,  (b)  methane,  (c)  carbon  monoxide  and  (d)  carbon  dioxide,  at  WGFR  of  9  and 
1  atm. 


hand,  when  hydrogen  is  removed  from  the  reactor,  it  does  not  only 
directly  affect  the  WGS  reaction  but  also  the  methanation  reaction, 
thus  inhibiting  the  last  reaction  to  a  higher  extent,  which  converts 
hydrogen  into  methane.  Therefore,  higher  hydrogen  yields  are  ob¬ 
tained  in  the  MR  than  in  the  SER  at  the  expense  of  higher  methane 
production  inhibition  (Figs.  3(b)  and  4(b)). 

As  already  stated,  when  there  is  no  methane  generation,  higher 
hydrogen  yields  are  reached,  being  that  the  temperature  of 
maximum  hydrogen  yield  depends  on  how  fast  methane  produc¬ 
tion  is  completely  inhibited.  In  the  SER,  this  maximum  occurs  at 
approximately  950  K  for  all  carbon  dioxide  removal  fractions  since 
carbon  dioxide  does  not  directly  interfere  in  the  methanation  re¬ 
action.  On  the  other  hand,  in  the  MR  the  temperature  at  which  the 
maximum  hydrogen  yield  is  obtained  decreases  for  increasing 
hydrogen  removal  fractions  (e.g.  900  K  for  40%  hydrogen  removal 
vs.  850  K  for  80%  hydrogen  removal).  This  emphasizes  once  again 
the  effect  that  hydrogen  removal  has  on  avoiding  more  intensively 
methane  formation  and  thus,  hydrogen  consumption. 

A  similar  behavior  is  observed  for  the  carbon  dioxide  yield  in 
both  reactors  (Figs.  3(d)  and  4(d)).  However,  on  the  contrary  to 
hydrogen,  the  yield  of  carbon  dioxide  starts  declining  at  a  tem¬ 
perature  bellow  the  one  that  completely  inhibits  methane  pro¬ 
duction,  for  the  reason  already  mentioned  in  the  previous 
section.  When  methane  production  is  completely  inhibited,  car¬ 
bon  dioxide,  in  the  same  way  as  hydrogen,  is  consumed  through 
RWGS.  Once  again,  in  a  SER  this  maximum  happens  approxi¬ 
mately  at  the  same  temperature  for  all  carbon  dioxide  removal 


fractions  while  in  the  MR  the  maximum's  temperature  decreases 
with  increasing  hydrogen  removal  fractions,  for  the  same  reasons 
presented  for  hydrogen.  Regarding  carbon  monoxide,  a  similar 
behavior  is  observed  in  both  SER  and  MR  (Figs.  3(c)  and  4(c)). 
Since  carbon  monoxide  is  consumed  through  WGS  and  metha¬ 
nation,  which  are  both  exothermic  reactions,  higher  tempera¬ 
tures  shift  the  opposite  reactions,  leading  to  higher  carbon 
monoxide  yields.  Finally,  when  higher  carbon  dioxide  or 
hydrogen  removals  are  applied,  more  carbon  monoxide  is 
consumed  through  WGS  and  so  lower  carbon  monoxide  yields  are 
obtained,  as  intended. 

3.2.  GSR  with  in  situ  hydrogen  and  carbon  dioxide  simultaneous 
removal 

In  this  section  the  effect  of  the  simultaneous  removal  of 
hydrogen  and  carbon  dioxide  during  GSR,  as  well  as  the  influence  of 
temperature,  WGFR  and  pressure  on  the  yield  of  hydrogen, 
methane,  carbon  monoxide  and  carbon  dioxide  is  accessed.  More¬ 
over,  an  analysis  of  the  boundaries  of  operating  conditions  for  solid 
carbon  formation  is  included  as  well.  During  the  simulations  for  the 
SEMR,  the  carbon  dioxide  removal  fraction  was  kept  at  0.99 
(approximate  removal  normally  obtained  before  the  breakthrough 
in  a  packed  bed)  while  the  hydrogen  removal  was  varied  between 
0  (SER  only)  and  0.99.  Once  again,  complete  glycerol  conversion 
and  no  coke  formation  were  observed  for  all  the  range  of  conditions 
analyzed. 
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Fig.  5.  Effect  of  temperature  and  hydrogen  removal  fraction  on  the  yield  of  (a)  hydrogen  and  (b)  methane  at  WGFR  of  9, 1  atm  and  carbon  dioxide  removal  fraction  of  0.99  in  a  SEMR. 


3.2 A.  Effect  of  temperature 

The  variation  of  the  yield  of  hydrogen  and  methane  with  tem¬ 
perature  and  hydrogen  removal  fraction  is  presented  in  Fig.  5.  Since 
the  variation  of  the  carbon  dioxide  yield  with  temperature  follows  a 
trend  very  similar  to  that  of  hydrogen,  its  graphic  is  not  presented 
here  for  brevity  reasons.  The  carbon  monoxide  yield  is  not  included 
here  as  well  since  only  residual  amounts  were  observed.  Constant 
WGFR  and  pressure  of  9  and  1  atm,  respectively,  were  used  since 
such  conditions  have  been  reported  in  many  works  13-16]  as  the 
most  favorable  conditions  for  GSR.  However,  their  effects  will  be 
analyzed  separately  in  a  later  section. 

As  can  be  seen  in  Fig.  5,  higher  temperatures  highly  enhance  the 
hydrogen  yield  and  decrease  the  methane  formation  because  of  the 
exothermic  nature  of  the  main  secondary  reaction,  methanation 
(Eq.  (4)),  which  takes  place  mainly  at  lower  temperatures. 
Regarding  the  hydrogen  removal  fraction,  it  clearly  enhances 
hydrogen  production,  especially  at  lower  temperatures.  For 
hydrogen  and  carbon  dioxide  removals  of  99%,  the  theoretical 
hydrogen  yield  of  7  is  obtained  for  all  the  range  of  temperatures 
tested.  However,  such  hydrogen  removal  in  a  membrane  reactor 
operating  at  1  atm  is  practically  impossible.  Therefore,  a  hydrogen 
removal  of  80%  [17  is  considered  to  be  much  more  realistic  while 
operating  at  1  atm.  Under  such  conditions  a  hydrogen  production  of 
5  and  7  mol/per  mole  of  glycerol  converted  at  600  I<  and  800  K, 
respectively,  was  obtained.  Compared  to  the  other  reactor  config¬ 
urations  (SER  and  MR),  the  SEMR  allows  higher  hydrogen  genera¬ 
tion,  especially  at  lower  temperatures.  Therefore,  from  the 


thermodynamic  point  of  view,  the  use  of  such  reactor  configuration 
would  be  much  more  beneficial  if  temperatures  around  600  I<  are 
used.  However,  this  decreases  reaction  kinetics  and  hydrogen 
permeation,  so  that  longer  reactors  and  thus,  longer  and/or  more 
permeable  membranes,  would  be  required. 

In  terms  of  methane  formation,  the  maximum  methane  pro¬ 
duction  of  0.5  moles/mole  of  glycerol  consumed  for  fCo2  =  0.99  and 
/h2  =  0.80  was  observed  at  600  K  (no  considerable  methane  at 
750  I<)  -  Fig.  5,  while  for  the  SER  with  fC02  =  0.99  and  for  the  MR 
with/n2  =  0.80,  maximum  methane  productions  of  1.29  (Fig.  3b) 
and  1.09  (Fig.  4b)  moles/mole  of  glycerol  converted  were  obtained, 
respectively,  for  the  same  temperature.  Regarding  carbon  dioxide, 
its  yield  varies  with  temperature  and  hydrogen  removal  fraction  in 
a  similar  way  as  hydrogen  (data  not  shown)  since  it  is  a  product  of 
the  same  reaction  that  origins  hydrogen  (WGS  (Eq.  (3))).  Moreover, 
since  the  base  case  here  already  considers  a  carbon  dioxide  removal 
fraction  as  high  as  0.99,  which  already  leads  to  a  very  similar 
behavior  between  hydrogen  and  carbon  dioxide  yields  in  the  SER 
(the  effect  of  carbon  dioxide  removal  outweighs  the  effect  of  tem¬ 
perature  -  Fig.  3),  then  it  is  expected  that  by  adding  hydrogen 
removal,  the  similarities  between  the  behavior  of  both  species 
yields  increase.  Finally,  it  is  again  worth  stressing  that  only  residual 
amounts  of  carbon  monoxide  (<0.01  mol/mol  of  converted  glyc¬ 
erol)  were  observed  during  the  simulations  for  all  the  range  of 
temperatures  and  hydrogen  removal  fractions. 

It  was  observed,  for  all  the  range  of  temperatures  and  hydrogen 
removal  fractions  under  which  methane  formation  occurs  that  the 


Fig.  6.  Effect  of  WGFR  and  hydrogen  removal  fraction  on  the  yield  of  (a)  hydrogen  and  (b)  methane  at  700  K,  1  atm  and  carbon  dioxide  removal  fraction  of  0.99  in  a  SEMR. 
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Fig.  7.  Effect  of  pressure  and  hydrogen  removal  fraction  on  the  yield  of  (a)  hydrogen  and  (b)  methane  at  700  K,  WGFR  of  9  and  carbon  dioxide  removal  fraction  of  0.99  in  a  SEMR. 


variation  of  temperature  or  hydrogen  removal  fraction  resulted  in 
variations  of  the  moles  of  hydrogen  4  times  higher  than  the  vari¬ 
ations  of  the  moles  of  methane  and  carbon  dioxide  (absolute  value). 
This  clearly  indicates  that  both  WGS  (Eq.  (3))  and  SRM  (reverse  of 
Eq.  (4))  control  the  GSR  process  in  the  SEMR.  On  the  other  hand, 
under  conditions  that  do  not  allow  methane  formation  (high 
hydrogen  removal  fractions  and  high  temperatures),  the  process  is 
solely  controlled  by  the  WGS  reaction. 

By  comparing  the  SEMR  with  the  other  three  reactor  configura¬ 
tions  in  terms  of  hydrogen  yield,  it  is  clear  that  the  SEMR  presents 
better  performance  than  both  the  SER  and  MR  and  much  better 
performance  than  the  TR.  In  fact,  while  in  the  TR,  at  700  K,  1  atm  and 
WGFR  of  9,  a  hydrogen  yield  of  2.18  was  attained,  in  the  SER  with 
fc o2  =  0.99,  MR  with  /Hz  =  0.80  and  SEMR  with  fCo2  =  0.99  and 
/h2  =  0.80,  hydrogen  yields  of  4.69,  5.66  and  6.92  were  observed, 
respectively,  under  the  same  conditions.  This  corresponds  to  an 
enhancement  of  the  hydrogen  yield  of  approximately  115%  when  the 
SER  with  fco2  =  0-99  is  used  instead  of  the  TR,  an  enhancement  of 
approximately  21%  when  the  MR  with  fHl  =  0.80  is  used  instead  of 
the  SER  with  fCo2  =  0.99  (enhancement  of  approximately  160% 
comparatively  to  the  TR)  and  an  increase  of  around  22%  when  the 
SEMR  with/co2  =  0.99  and  fHl  =  0.80  is  used  instead  of  the  MR  with 
/h2  =  0.80  (enhancement  of  217%  and  47%  comparatively  to  the  TR 
and  the  SER  with  fCo2  =  0.99,  respectively). 

3.2.2.  Effect  of  water-to-glycerol  feed  ratio 

The  effect  of  both  hydrogen  removal  fraction  and  WGFR  on  the 
yield  of  hydrogen  and  methane  was  target  of  investigation  as  well. 
Such  study  is  depicted  in  Fig.  6.  The  WGFR  was  varied  between  3 
(stoichiometric  value)  and  9,  which  has  been  considered  by  many 
authors  as  the  most  favorable  value,  as  already  mentioned.  Higher 
values  were  not  considered  because  this  would  imply,  in  practice, 
excessive  costs  for  water  vaporization. 

In  terms  of  hydrogen  production,  the  increase  of  the  WGFR 
enhances  the  production  of  hydrogen,  especially  at  lower  hydrogen 
removal  fractions.  By  increasing  the  water  content  in  the  reaction 
environment,  the  thermodynamic  equilibrium  of  both  GSR  and 
SRM  (reverse  of  Eq.  (4))  is  dislocated,  according  to  the  Le  Chatelier's 
principle,  towards  the  conversion  of  the  excess  of  water  into 
hydrogen.  Thus,  the  methane  yield  decreases,  as  can  be  seen  in 
Fig.  6(b).  Regarding  the  carbon  dioxide  yield,  its  behavior  is  very 
similar  to  hydrogen  (data  not  shown)  because  of  the  reasons 
already  stated  in  the  previous  section.  Once  again,  only  residual 
amounts  of  carbon  monoxide  were  observed.  The  process  is 
controlled  by  the  reactions  in  Eqs.  (2)— (4)  for  all  the  WGFRs 
involved.  When  there  is  no  methane,  only  Eqs.  (2)  and  (3)  are 
relevant. 


Concerning  the  effect  of  increasing  hydrogen  removal  fractions, 
it  is  clear  that  higher  removal  fractions  of  hydrogen  lead  to  higher 
hydrogen  production  and  methane  production  inhibition.  More¬ 
over,  this  enhancement  is  more  pronounced  for  less  favorable 
conditions  (WGFR  =  3). 

In  sum,  considering  again  the  case  in  which  fCo2  =  0.99  and 
/h2  =  0.80,  when  the  WGFR  increased  from  3  to  9,  the  hydrogen 
yield  increased  from  4.57  to  6.92  moles/mole  of  converted  glycerol, 
while  the  yield  of  methane  decreased  from  0.61  to  0.02. 

If  higher  WGFRs  had  been  tested,  higher  hydrogen  yields  would 
have  been  obtained,  especially  at  lower  hydrogen  removal  frac¬ 
tions.  However,  such  water  contents  are  considered  to  be  disad¬ 
vantageous  due  to  the  high  vaporization  costs  involved  at  industrial 
scale.  Moreover,  even  if  they  allow  avoiding  coke  formation,  high 
water  contents  may  lead  to  catalyst  and  sorbent  deactivation  due  to 
pore  blocking  [18  .  On  the  contrary,  lower  WGFRs  are  not  so 
problematic  from  the  economic  point  of  view  since  the  glycerol 
produced  through  the  biodiesel  production  process  is  composed  by 
low  amounts  of  water  (3.2%  water,  40%  fatty  matter,  33%  glycerol, 
23%  methanol  and  3.8%  ash)  19].  Therefore,  if  lower  WGFRs  are 
employed,  less  water  has  to  be  added.  Even  if  low  WGFRs  are  not  so 
advantageous  in  terms  of  hydrogen  production,  the  use  of  a  hybrid 
reaction-separation  process  like  the  one  here  described  allows 
reducing  the  constraints  imposed  by  such  low  water  contents. 

3.2.3.  Effect  of  pressure 

Finally,  the  effect  of  pressure  and  hydrogen  removal  fraction  on 
the  yield  of  hydrogen  and  methane  was  analyzed,  as  can  be  seen  in 
Fig.  7.  It  can  be  observed  that  when  higher  pressures  are  used,  the 
yield  of  hydrogen  suffers  a  negative  effect,  except  for  very  high 
hydrogen  removal  fractions  where  constant  hydrogen  yields  are 
observed.  Such  effect  is  once  again  related  to  the  Le  Chatelier's 
principle,  because  during  GSR  the  total  number  of  moles  varies  (Eq. 
(1)).  Therefore,  when  higher  pressures  are  applied,  the  thermody¬ 
namic  equilibrium  will  be  shifted  towards  the  production  of  a  lower 
total  number  of  moles.  However,  once  the  glycerol  conversion  was 
complete  for  all  the  range  of  pressures  used,  the  decrease  of  the 
hydrogen  yield  cannot  be  associated  to  GSR.  Instead,  it  is  associated 
to  the  enhancement  of  the  methanation  reaction  (Eq.  (4)),  which 
consumes  hydrogen  and  yields  methane.  As  already  stated,  this  is 
not  valid  for  very  high  hydrogen  removal  fractions,  where  there  is 
no  methane.  While  at  1  atm,  and  considering  fCo2  =  0.99  and/n2  = 
0.80,  the  hydrogen  and  methane  yields  were  6.92  and  0.02, 
respectively,  at  5  atm  the  hydrogen  yield  decreased  to  6.06  and  the 
methane  one  increased  to  0.23. 

In  terms  of  carbon  dioxide,  its  yield  varies  with  pressure  in  a  very 
similar  way  as  hydrogen  does  (data  not  shown),  for  the  reasons 
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already  mentioned.  The  production  of  carbon  monoxide  was  once 
again  negligible.  Also,  both  RWGS  (reverse  of  Eq.  (3))  and  metha- 
nation  (Eq.  (4))  control  the  process  for  all  the  range  of  pressures 
analyzed,  except  when  the  formation  of  methane  is  completely 
inhibited.  In  this  last  case,  WGS  is  the  only  relevant  reaction. 

Regarding  the  effect  of  the  hydrogen  removal  fraction,  analo¬ 
gously  to  the  previous  cases,  the  hydrogen  yield  is  favored  by 
increasing  hydrogen  removals,  while  methane  suffers  an  opposite 
effect.  The  higher  enhancement  of  the  hydrogen  yield  is  observed  at 
the  less  thermodynamically  favorable  pressure  (5  atm).  It  is, 
nonetheless,  important  to  mention  that  the  approach  used  here 
ignores  the  effect  of  pressure  on  the  removal  of  both  carbon  dioxide 
and  hydrogen,  as  well  as  in  the  reaction  kinetics.  It  is  known  that 
pressure  influences  the  driving  force  for  carbon  dioxide  sorption 
and  for  hydrogen  permeation  through  a  hydrogen  perm-selective 
membrane.  Thus,  higher  pressures  lead  to  higher  carbon  dioxide 
and  hydrogen  removal  fractions,  meaning  that  if  high  pressures  are 
used  in  the  SEMR,  high  hydrogen  removals  are  attained  and 
inherently  higher  hydrogen  yields.  Moreover,  in  the  case  of 
fc o2  =  0.99  and/n2  =  0.80,  when  the  pressure  increases  from  1  to 
5  atm,  the  hydrogen  yield  at  equilibrium  decreases  much  less  than 
for  lower  hydrogen  removals,  namely  for  fHl  =  0  (SER).  Conse¬ 
quently,  the  use  of  high  pressures  becomes  advantageous  after  all 
due  to  the  membrane  effect. 

3.2.4.  Coke  formation 

The  formation  of  coke  was  target  of  investigation  as  well,  since  it 
is  an  undesired  product  that  may  affect  the  performance  of  the 
catalyst,  sorbent  and  membrane  during  actual  operation  of  the 
SEMR.  It  was  observed  that  the  formation  of  coke  in  the  TR  in¬ 
creases  with  the  decrease  of  temperature  (and  so  a  higher  WGFR 
would  be  required  to  completely  avoid  coke  formation)  during  GSR, 
however  at  600  I<  and  1  atm  a  WGFR  of  4  is  high  enough  to 
completely  avoid  coke  formation.  When  carbon  dioxide  or 
hydrogen  were  removed  from  the  reaction  medium,  opposite  ef¬ 
fects  were  observed.  While  the  formation  of  coke  decreases 
comparatively  to  the  TR  when  carbon  dioxide  is  removed  from  the 
reaction  medium,  when  hydrogen  is  removed  more  coke  is  pro¬ 
duced.  For  a  carbon  dioxide  removal  of  99%,  a  WGFR  of  1  is  enough 
to  completely  avoid  coke  formation  at  temperatures  between  600 
and  800  K  and  1  atm.  On  the  other  hand,  when  99%  of  hydrogen  is 
removed,  a  WGFR  of  4.5  is  necessary  to  avoid  carbon  formation  in 
the  same  temperature  range  and  atmospheric  pressure.  Therefore, 
it  allows  us  to  conclude  that  most  of  coke  is  formed  from  methane 
through  the  reaction  in  Eq.  (8)  (the  low  WGFRs  allow  the  formation 
of  methane,  even  at  hydrogen  and  carbon  dioxide  removals  of  99%), 
being  that  the  reactions  in  Eqs.  (7),  (9)  and  (10)  are  not  compatible 
with  the  observed  behavior  and  so  they  can  be  neglected. 

It  is  therefore  expected  that  when  both  products  are  removed 
from  the  reaction  medium  in  the  SEMR,  an  equilibrium  between 
both  opposite  effects  is  reached.  That  was  in  fact  observed,  being  that 
at  /co2  =  0.99  and  fHl  =  0.80  a  WGFR  of  2.5  was  enough  to 
completely  inhibit  coke  formation  for  all  the  range  of  temperatures 
and  pressures  tested.  Even  for  a  removal  of  99%  of  both  products,  the 
stoichiometric  WGFR  of  3  was  sufficient  to  avoid  coke  formation  in 
all  the  temperature  and  pressure  range.  Therefore,  it  would  be 
possible,  from  the  thermodynamic  point  of  view,  to  operate  a  SEMR 
for  GSR  under  the  temperature,  WGFR  and  pressure  ranges  analyzed 
in  this  study  without  having  any  problems  with  coke  formation. 

4.  Conclusions 

A  thermodynamic  analysis  was  carried  out  on  the  GSR  for 
hydrogen  production  in  a  SEMR  with  in  situ  hydrogen  and  carbon 
dioxide  simultaneous  removal.  At  800  K,  WGFR  of  9,  1  atm  and 


fco2  =  0.99  and  fa2  =  0.80,  the  hydrogen  yield  was  7  (maximum 
possible  value).  Regarding  methane  and  carbon  monoxide,  no 
considerable  production  was  observed  under  such  conditions.  Even 
though  the  SEMR  presents  higher  hydrogen  production  than  the 
SER  and  the  MR  and  much  higher  hydrogen  production  than  the  TR, 
the  difference  is  not  very  significative  at  higher  temperatures 
compared  to  the  two  first  configurations.  However,  this  difference 
becomes  much  more  meaningful  with  the  decrease  of  temperature. 
At  600  K,  a  hydrogen  production  of  5  moles/mole  of  converted 
glycerol  was  observed  in  the  SEMR  with  fCo2  =  0.99  and  =  0.80, 
while  hydrogen  productions  of  1.83  and  2.62  moles/mole  of  con¬ 
verted  glycerol  were  obtained  in  the  SER  and  MR  with  the  same 
removals,  respectively.  Therefore  the  advantages  of  the  SEMR  are 
more  noticeable  at  lower  temperatures  and  so,  this  new  hybrid 
reactor  configuration  can  be  seen  as  a  possible  solution  to  lower  the 
high  operating  temperatures  necessary  for  the  endothermic  GSR. 
However,  a  proper  judgment  has  to  be  done,  because  at  lower 
temperatures  longer  reactors  and  longer  and/or  more  permeable 
membranes  would  be  necessary  to  compensate  the  decreased  re¬ 
action  kinetics  and  permeation  rate. 

In  terms  of  WGFR  and  pressure  effects,  it  is  clear,  from  ther¬ 
modynamic  point  of  view,  that  higher  WGFRs  and  lower  pressures 
favor  the  production  of  hydrogen  in  the  SEMR.  Finally,  it  was 
observed  that  no  coke  was  formed  at  WGFRs  higher  than  3. 

Acknowledgments 

Joel  M.  Silva  acknowledges  the  Doctoral  Program  in  Chemical 
and  Biological  Engineering  from  FEUP  and  the  FCT  for  the  PhD  grant 
(PD/BD/52625/2014).  M.A.  Soria  is  grateful  to  the  FCT  for  the 
postdoctoral  grant  (SFRH/BPD/88444/2012). 

Appendix  A.  Supplementary  data 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
dx.doi.org/10.1016/j.jpowsour.2014.09.093. 

References 

[1]  H.  Chen,  T.  Zhang,  B.  Dou,  V.  Dupont,  P.  Williams,  M.  Ghadiri,  Y.  Ding,  Int.  J. 
Hydrogen  Energy  34  (2009)  7208—7222. 

[2]  Y.  Li,  W.  Wang,  B.  Chen,  Y.  Cao,  Int.  J.  Hydrogen  Energy  35  (2010)  7768-7777. 

[3]  X.  Wang,  N.  Wang,  M.  Li,  S.  Li,  S.  Wang,  X.  Ma,  Int.  J.  Hydrogen  Energy  35 
(2010)  10252-10256. 

[4]  C.V.  Miguel,  A.  Mendes,  S.  Tosti,  L.M.  Madeira,  Int.  J.  Hydrogen  Energy  37 
(2012) 12680-12687. 

[5]  J.M.  Silva,  M.A.  Soria,  L.M.  Madeira,  submitted  for  publication. 

[6]  A.C.C.  Chang,  W.-H.  Lin,  K.-H.  Lin,  C.-H.  Hsiao,  H.-H.  Chen,  H.-F.  Chang,  Int.  J. 
Hydrogen  Energy  37  (2012)  13110—13117. 

[7]  K.-H.  Lin,  W.-H.  Lin,  C.-H.  Hsiao,  H.-F.  Chang,  A.C.C.  Chang,  Int.  J.  Hydrogen 
Energy  37  (2012)  13770-13776. 

[8]  A.  Hanif,  S.  Dasgupta,  S.  Divekar,  A.  Arya,  M.O.  Garg,  A.  Nanoti,  Chem.  Eng.  J. 
236  (2014)  91-99. 

[9]  C.V.  Miguel,  R.  Trujillano,  V.  Rives,  M.A.  Vicente,  A.F.P.  Ferreira,  A.E.  Rodrigues, 
A.  Mendes,  L.M.  Madeira,  Sep.  Purif.  Technol.  127  (2014)  202—211. 

[10]  P.N.  Sutar,  P.D.  Vaidya,  A.E.  Rodrigues,  Chem.  Eng.  Technol.  33  (2010) 
1645-1649. 

[11]  C.K.  Cheng,  S.Y.  Foo,  A.A.  Adesina,  Ind.  Eng.  Chem.  Res.  49  (2010) 
10804-10817. 

[12]  S.  Adhikari,  S.  Fernando,  S.R.  Gwaltney,  S.D.  Filip  To,  R.  Mark  Bricka, 
P.H.  Steele,  A.  Haryanto,  Int.  J.  Hydrogen  Energy  32  (2007)  2875—2880. 

[13]  X.  Wang,  S.  Li,  H.  Wang,  B.  Liu,  X.  Ma,  Energy  Fuels  22  (2008)  4285-4291. 

[14]  H.  Chen,  Y.  Ding,  N.T.  Cong,  B.  Dou,  V.  Dupont,  M.  Ghadiri,  P.T.  Williams, 
Renew.  Energy  36  (2011)  779—788. 

[15]  S.  Adhikari,  S.  Fernando,  A.  Haryanto,  Energy  Fuels  21  (2007)  2306—2310. 

[16]  N.  Hajjaji,  A.  Chahbani,  Z.  Khila,  M.-N.  Pons,  Energy  64  (2014)  473—483. 

[17]  D.  Mendes,  S.  Sa,  S.  Tosti,  J.M.  Sousa,  L.M.  Madeira,  A.  Mendes,  Chem.  Eng.  Sci. 
66  (2011)  2356-2367. 

[18]  Z.  Martunus,  A.D.  Helwani,  J.  Wiheeb,  M.R.  Kim,  Othman,  Int.  J.  Greenhouse 
Gas  Control  7  (2012)  127-136. 

[19]  M.  Slinn,  K.  Kendall,  C.  Mallon,  J.  Andrews,  Bioresour.  Technol.  99  (2008) 
5851-5858. 


